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ABSTRACT 

The non thermal radiation observed from a handful of clusters of galaxies (CG) 
is the proof that particle acceleration occurs in the intra cluster medium (ICM). 
It is often believed that shock surfaces associated with either mergers of CG, or 
with the cosmological inflow of matter onto clusters during structure formation 
may be the sites for acceleration. We discuss here the effectiveness of shock 
acceleration in the ICM, stressing that merger related shocks are typically 
weak, at least for the so-called major mergers. We investigate the implica- 
tions of shock strengths for gamma ray emission from single CG and for their 
detectability with future gamma ray satellites (such as GLAST and AGILE) 
and ground based Cherenkov telescopes. We also discuss the contribution of 
clusters to the extragalactic diffuse gamma ray background (EDGRB). 



1 Introduction 



Clusters of galaxies are powerful X-ray sources, with luminosities in the range 
10^'^ — W^^erg s^^ . This thermal radiation is emitted by an hot (~ 10* if), 
tenuous 10~^cm^^), highly ionized intra cluster gas that fills all the cluster 
volume and accounts for a dominant fraction of the cluster baryons D . Com- 
pelling evidence for the presence of relativistic particles and magnetic field in 
the intra cluster medium (hereafter ICM) is given by the observation of a diffuse 
synchrotron radio emission in a consistent fraction of the rich clusters ^ . This 
gives us a direct proof of the presence of GeV electrons and ^ fiG magnetic 
field. The same electrons responsible for the radio emission can upscatter the 
cosmic microwave background (CMB) photons in the X-ray band via inverse 
Compton scattering (ICS) providing an explanation for the hard X-ray excess 
observed in a few clusters ^ ^ . Different mechanisms, such as non thermal 
Bremsstrahlung from a population of supra thermal electrons, have also been 
proposed 0*. 

Clusters are also expected to be potential gamma ray sources, the main 
emission processes being inverse Compton scattering off the cosmic microwave 
background photons (CMB) from ultra relativistic electrons and the decay of 
neutral pions produced in the interactions between cosmic ray protons and 
ICM protons ^ ' . In this paper we make predictions on the number of clusters 
that future gamma ray telescopes (both from space and earth) will be able to 
detect and on the contribution of clusters to the extragalactic diffuse gamma 
ray background. We focus our attention on the photons produced by electrons 
accelerated at shocks that develop during the process of large scale structure 
formation. Protons are also accelerated at shocks, they cannot escape the clus- 
ter volume £3 and they are re-energized each time a new merger event occurs 
IS). It follows that the spectrum and the energy density of protons depend 
on the merger history of the cluster and an hadronic feature might appear on 
top of the gamma ray emission generated by inverse Compton scattering of 
electrons . 

2 Shock acceleration during structure formation 

The ICM is heated to the observed high temperature by shock waves that 
naturally form during the process of large scale structure formation '^^^ . Dif- 



fusive acceleration can take place at these shocks so that a small fraction of 
the particles crossing the shock discontinuity can be energized up to relativistic 
energies. The maximum energy for the accelerated particles can be evaluated 
equating the acceleration time with the energy loss time. Under the reasonable 
assumption that during the acceleration process the residence time of the par- 
ticle in the upstream region is longer than the residence time downstream, and 
neglecting a weak dependence on the shock Mach number we can approximate 
the acceleration time as ; 

where D and u are the upstream diffusion coefficient and the upstream fluid 
velocity. Under these hypothesis the acceleration time depends only on up- 
stream quantities, unaffected by the presence of the shock. Assuming a Bohm 
diffusion coefficient we obtain: 

Tacc{E) - O.W{iiG)-'EiGeV)i—^ -) years. (2) 

f 0*cms~^ 

For electrons the relevant loss channel is ICS on the CMB radiation character- 
ized by a time scale: 
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l(fE{GeV)-^{l + z)-^years (3) 



where Ucaib is the CMB energy density and 7 is the electron Lorentz factor. 
Equating 121 and 13 we obtain the maximum energy for the electrons: 

E„ra. ~ 57B{iiG)^/^i—^ -) TeV (4) 

corresponding to a maximum energy for the upscattered CMB photons equal 
to: 

4 11 
i?— = -^hcMB - 7.5S(MG)(-g TeV (5) 

which falls in the very high energy gamma ray band. It is important to stress, 
however, that the Bohm diffusion coefficient is the smallest possible one, and 
that a different choice of the diffusion coefficient could result in a too low max- 
imum energy of the electrons, so that the electron population would be unable 
to upscatter the CMB photons up to gamma ray energies. In the following we 
use for the maximum energy the value obtained in eq. [S] 



The spectra of particles accelerated at shocks are power laws in momen- 
tum of the form N{p) oc where the slope a is comrected to the shock Mach 
number M through the well known relation: 

3 Large scale structure formation and associated shocks. 

According to the standard theory of large scale structure formation clusters 
form hierarchically through mergers of smaller subunits. Semi-analytical mod- 
els describing this process can be found in the literature _ This for- 
malism provides us the comoving density of objects with given mass M at a 
redshift z, n{M, z), and the rate at which clusters of mass Mi and M2 merge 
together at a fixed redshift z, R{Mi, M2, z). Sampling these distributions it 
is possible to simulate the whole merger history of a cluster with a certain 
present mass. Each time a merger occurs we can estimate the relative velocity 
Vr between the two involved clusters using energy conservation EI IS) ; 

_ GM1M2 ^ 1^^^2 ^ GM1M2 
Rl+ R2 2 2i?(a 

where Mi and Ri are the masses and virial radii of the two colliding clusters, 
Mr is the reduced mass and Rta is the turnaround radius of the system. The 
shock Mach number is the ratio between Vr and the sound speed in the ICM, 
given by the virial theorem: 
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= 7(7-1)^- (8) 



In figure (left panel) the Mach numbers M of the shocks are plotted as a 
function of the ratio between the masses of the two merging clusters. This 
result has been obtained simulating 500 different realizations of a merger tree 
for a cluster with present mass 10^^ solar masses In the right panel the 
cumulative distribution of Mach numbers is shown. 

Fig. ntells us that merger shocks are often weak, following a distribution 
strongly peaked at '--^ 1.4 and consequently are unable to accelerate rela- 
tivistic particles in an efficient way (accelerated spectra are very steep power 
law with a ~ 6). Only ^ 6% of the shocks have Mach number greater than 3 
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Figure 1: Merger shock Mach numbers as a function of the ratio between the 
masses of the two colHding chisters (left panel) and Mach number cumulative 
distribution (right panel) . The upper strip in the left panel refers to the shocks 
that propagate in the smaller merging cluster while the lower strip refers to the 
shocks propagating in the bigger, hotter cluster. 



(corresponding to a > 2.4), necessary to be consistent with the spectral shape 
of cluster radio halos, and these strong shocks form only during minor mergers 
(collisions between clusters with very different masses). On the contrary radio 
halos are found in a consistent fraction (~ 40%) of massive clusters and there 
seems to be a correlation between the presence of the halo and the signature 
of a recent major merger EJ . We can conclude that radio halos are unlikely to 
be produced by electrons directly accelerated at merger shocks ^ . 

Besides merger shocks, that develop in the ICM of virialized objects, we 
consider here also accretion shocks that form due to the secondary infall of non 
virialized matter onto an already formed cluster . These shocks propagate 
in a cold medium and are strong by definition, so that the spectra of the 
accelerated particles are fixed to N{p) oc p~^. We assume here that a spherical 
shock forms at the virial radius of each cluster and that the cold infalling 
matter (whose density is taken here equal to the average baryon density of the 
universe) flows at the free fall velocity. 



4 Source counts and the contribution to the gamma ray background 



Recently an association between clusters of galaxies and unidentified EGRET 
sources was proposed by different authors Uiil I1Z|I. The physical plausibility 
and the statistical significance of such an association were strongly questioned 
in "iSi) ^ and the lack of association was also found in . Future gamma ray 
space telescopes, such as AGILE and GLAST, will observe the whole sky with 
a better sensitivity compared to EGRET. We make here predictions on the 
number of objects that these telescopes will be able to detect. 

The number of accreting clusters with gamma ray flux greater than a 
fixed value can be evaluated as follows: 



r°° dV 

N{>Fu,n)= I dz— dMn{M,z) (9) 



dz 



where dV is the comoving volume between z and z + dz, n{M, z) is given in , 
and M(Fiirm z) is the mass of a cluster accreting at redshift z whose flux is Fum- 
For merging clusters a similar, slightly more complicated expression, which 
includes also the merger rate i?(Afi, M2, z), holds. Gamma ray luminosities 
are evaluated assuming that a small fraction 77 ^ 5% of the kinetic energy 
flowing across the shock is converted in relativistic electrons, which upscatter 
the CMB photons in the gamma ray energy band. 

The results, as obtained in are shown in fig |2 (left panel) and can be 
summarized in this way: AGILE and GLAST will be able to observe ^ 10 and 
a few tens of clusters respectively (both merging and accreting), while no one 
of the unidentified EGRET sources is expected to be associated with a cluster, 
in agreement with li^. 

The superposition of gamma ray emission from single clusters results in a 
diffuse background. In was firt proposed that the CMB photons, scattered 
up to gamma ray energies by relativistic electrons accelerated at the shocks that 
form during large scale structure formation, could account for the whole ex- 
tragalactic diffuse gamma ray back grou nd, as estimated sub tracting a galactic 
emission model to the EGET data l^iJ . We evaluate in this contribution 
using the formalism described in the previous sections. Our results are shown 
in fig. 121 (right panel). The shaded region corresponds to the EGRET data 
and the solid line that roughly fits the data is the prediction made in 129. 
Our estimates of the contribution to the background from accreting and merg- 



Figure 2: Number of accreting (solid line) and merging (dashed line) clusters 
with gamma ray flux greater than F (left panel) . In the right panel are shown 
our predictions for the extragalactic diffuse gamma ray background (see text 
for details). 



ing clusters are represented by the dotted and dashed lines respectively. The 
dot-dashed line represents a more realistic estimate for the diffuse emission due 
to merging clusters, obtained considering only objects with mass greater than 
(the typical mass of a group of galaxies). According to our estimate, 
forming structures in the universe can account only for 10% of the extra- 
galactic background a nd this i s in agreement with recent results from detailed 
numerical simulations . 



5 Very high energy gamma ray emission 

The maximum energy of the upscattered CMB photons falls in the TeV range 
(see eqISJ so that we can wonder if clusters could be detected by ground based 
Cherenkov telescopes. In fig. Elwe plot the sensitivities for a generic telescope 
as calculated in for 100 hours of observation with an array consisting 

of 10 cells (pointlike and 1° sources) together with the expected spectra for 
three different cases (top to bottom): 1) a merger between two clusters with 
masses 10^^ and 10^'^ Mq, 2) a IO^^Mq accreting cluster with a magnetic field 
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Figure 3: Gamma ray emission in the 100 GeV - 10 TeV region. The thick 
sohd Hnes represent the sensitivities of a cherenkov array for point sources 
(lower curve) and extended sources (upper curve). The predicted gamma ray 
fluxes from a Coma-like cluster with and without absorption of the infrared 
background are plotted as dashed and solid lines respectively. 



in the upstream region equal to O.l/^G, 3) the same as 2) but with a ten times 
lower magnetic field. All the considered clusters are at a distance of lOOMpc. 
The thin continuous and dotted lines refer to the spectra calculated with and 
without taking into account the absorption on the infrared background 1201 . 

The apparent size of a 10^^ Mq accreting cluster at a distance of lOOMpc 
is comparable or even greater than the field of view of Cherenkov telescopes 
a few degrees), while the size of the emitting region for the merger considered 
in fig. |31should be roughly 1°, making it marginally detectable. 



6 Conclusions 



In Isji we first proposed that merger shocks are often weak and consequently 
unable to accelerate efficiently relativistic particles. We developed a recipe 
to evaluate the shock Mach numbers as a function of the masses of the two 
merging objects finding that shock acceleration can be effective only during 
minor mergers. We extended our analysis in where we included also the 
effects of strong accretion shocks. 

Ultrarelativistic electrons can be accelerated at such strong shocks and 
they can upscatter the CMB photons to very high gamma ray energies, making 
clusters of galaxies potential GeV and TeV gamma ray sources. We estimate 
that the future gamma ray telescopes AGILE and GLAST will observe ~ 10 and 
~ 50 clusters respectively ^ , while we do not expect any association between 
clusters and unidentified EGRET sources, in agreement with the observational 
results presented in We also estimated the contribution from large scale 

structure formation to the EDGRB to be at most equal to ~ 10% of the 
observed one. 

Only massive, nearby merging clusters could be detected by future Cherenkov 
arrays, while the detection of accreting clusters seems to be more problematic 
due to the very extended apparent size of these sources. 
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